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ABSTRACT 

We report extensive radio and X-ray observations of SN 2003bg whose spec¬ 
troscopic evolution shows a transition from a broad-lined Type Ic to a hydrogen- 
rich Type II and later to a typical hydrogen-poor Type Ibc. We show that the 
extraordinarily luminous radio emission is well described by a self-absorption 
dominated synchrotron spectrum while the observed X-ray emission at f ~ 30 
days is adequately £t by Inverse Compton scattering of the optical photons off 
of the synchrotron emitting electrons. Our radio model implies a sub-relativistic 
ejecta velocity, v ~ 0.24c, at fo ~ 10 days after the explosion which emphasizes 
that broad optical absorption lines do not imply relativistic ejecta. We hnd that 
the total energy of the radio emitting region evolves as i? ~ 7.3 x 
erg assuming equipartition of energy between relativistic electrons and magnetic 
helds (ce = eB = 0.1). The circumstellar density is well described by a stellar 
wind prohle with modest (factor of ~ 2) episodic density enhancements which 
produce abrupt achromatic flux variations. We estimate an average mass loss 
rate oi M ^ 3 x 10“^ Mq yr“^ (assuming a wind velocity of = 10^ km s“^) 
for the progenitor, consistent with the observed values for Galactic Wolf-Rayet 
stars. Comparison with other events reveals that ~ 50% of radio supernovae show 
similar short timescale flux variations attributable to circumstellar density irreg¬ 
ularities. Specifically, the radio light-curves of SN 2003bg are strikingly similar 
to those of the Type Ilb SN 2001ig, suggestive of a common progenitor evolu¬ 
tion for these two events. Based on the relative intensity of the inferred density 
enhancements, we conclude that the progenitors of SNe 2003bg and 2001ig expe¬ 
rienced quasi-periodic mass loss episodes just prior to the SN explosion. Finally, 
this study emphasizes that abrupt radio light-curve variations cannot be used 
as a reliable proxy for an engine-driven explosion, including off-axis gamma-ray 
bursts. 
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1. Introduction 

Accounting for ~ 10% of the nearby supernova population, Type Ibc supernovae (here¬ 
after SNe Ibc) are identihed by their lack of spectroscopic hydrogen and silicon features (see 
Filippenko 1997 for a review). Recognized as a rare subclass of core-collapse supernova ~ 20 
years ago (Elias et al. 1985), SNe Ibc have recently enjoyed a revitalized interest. Beginning 
with the discovery of Type Ic SN 1998bw in temporal and spatial coincidence with gamma- 
ray burst (GRB) 980425 (Galama et al. 1998; Pian et al. 2000), we now know that most 
long-duration GRBs (e.g. Stanek et al. 2003) and X-ray flashes (XRFs; Soderberg et al. 
2005b) are associated with SNe Ibc. However, radio observations of local SNe Ibc indicate 
that the inverse is not true; there is a strict limit of < 10% on the fraction of SNe Ibc that 
could be accompanied by a GRB or XRF (Soderberg et al. 2006). The lack of hydrogen 
features in both SNe Ibc and GRB/XRF-associated SNe imply that they represent the ex¬ 
plosion of massive stripped-core progenitors (Wheeler & Levreault 1985; Woosley, Eastman 
& Schmidt 1999). 

The popular models for SNe Ibc progenitors include massive Wolf-Rayet (WR) stars 
that eject their envelopes through strong dense winds, and close binary systems where the 
progenitor star is stripped of its hydrogen-rich layer by the companion (Ensman & Woosley 
1988). Despite dedicated archival searches, the progenitors of SNe Ibc are still poorly con¬ 
strained by pre-explosion images. The best photometric constraints are currently associated 
with SN 2002ap and SN 2004gt which exclude WR progenitors in the top ~ 30% and ~ 50% 
of the population, respectively (Smartt et al. 2002; Gal-Yam et al. 2005; Maund, Smartt & 
Schweizer 2005). While intriguing, these archival observations lack the sensitivity to clearly 
discriminate between the WR and binary progenitor models. 

By studying the circumstellar medium around the supernova, however, it is possible 
to place independent constraints on the progenitor. Observations show that Galactic Wolf- 
Rayet stars are embedded in wind stratihed media (Ghu 2002) while binary systems are 
associated with disrupted circumstellar media, possibly including an outflow during a com¬ 
mon envelope phase (Podsiadlowski, Joss & Hsu 1992). 

Radio observations provide the most direct probe of the circumstellar density structure 
around supernovae (Weiler et al. 1986; Ghevalier 1998). By observing the dynamical in¬ 
teraction of the ejecta with the surrounding medium we are able to map out the mass loss 
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history of the progenitor star. The identihcation of irregular density prohles may therefore 
distinguish between Wolf-Rayet and binary progenitor systems. 

Here we present extensive radio observations of SN 2003bg, discovered as part of our 
ongoing radio survey of local Type Ibc supernovae. The peculiar SN 2003bg spectroscopi¬ 
cally evolved from a broad-lined Type Ic to a hydrogen-rich Type II and later to a typical 
Type Ibc (Hamuy et ai, in prep), thereby bridging the hydrogen-rich and poor divisions 
of the core-collapse classification system. Our densely-sampled radio light-curves show that 
the extraordinarily luminous radio emission for a SN Ibc is characterized by episodic short- 
timescale variations. We show that these variations are well described by abrupt density 
enhancements in the circumstellar medium. Comparison with other radio supernovae shows 
that ~ 50% of all well-studied events similarly show evidence for abrupt light-curve varia¬ 
tions. We review the mass loss evolution observed (and inferred) for the progenitors of these 
core-collapse SNe. Using the observed radio properties for the peculiar SN 2003bg, we place 
constraints on the nature of its progenitor system. 

The organization of this paper is as follows: observations from the Very Large Array 
(VLA) and the Chandra X-ray Observatory (CXO) are described in §2. Preliminary estimates 
of the energy, velocity and density of the radio emitting region are presented in §3. Modeling 
of the radio light-curves are presented in §4 and §5 while a discussion of the radio polarization 
mechanism follows as §6. Our modeling of the X-ray emission is discussed in §7. In §8 we 
present a compilation of radio supernovae with abrupt light-curve variations and review 
their circumstellar irregularities. Finally, in §9 we discuss the possible causes of the density 
enhancements surrounding SN 2003bg and the implications for the progenitor system. 


2. Observations 

SN2003bg was optically discovered on 2003 February 25.7 UT, offset 16”.3 W and 
24”.6 S from the center of host galaxy, MCG -05-10-15, at d ^ 19.6 Mpc (Wood-Vasey & 
Chassagne 2003). Early spectroscopy on 2003 February 28 UT indicated that SN 2003bg 
was a peculiar Type Ic supernova with broad optical absorption lines, indicative of fast 
photospheric velocities and comparable to those seen in SN 1998bw (Filippenko & Chornock 
2003). Soon thereafter, SN 2003bg developed strong, broad Ha emission with a P-Cygni 
absorption component indicating an expansion velocity of —17,110 km s“^ (Hamuy, Phillips 
& Thomas-Osip 2003). P-Cygni components were additionally observed for H/5 and Hy, 
prompting the reclassihcation of SN 2003bg as an unusual Type H supernova. Nebular 
spectra taken several months later revealed that all evidence of hydrogen had disappeared, 
suggestive that the SN had transitioned back to a Type Ibc event (Hamuy et ah, in prep). 
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The peak time of the optical light-curves suggests that SN 2003bg was discovered within 
just a few days of explosion (Hamuy, private communication). Throughout this paper we 
therefore assume an approximate explosion date of 2003 February 22 UT. 


2.1. Very Large Array Data 

Motivated by the spectroscopic similarity of SN 2003bg to the radio luminous SN 1998bw, 
we initiated radio observations of SN 2003bg on 2003 March 4.07 UT with the Very Large 
Array^ (VLA). At 8.46 GHz we detected a radio source coincident with the optical position 
at a(J2000) = 04^10“59.42^ 5(J2000) = -31°24'50.3" (±0.1 arcsec in each coordinate) with 
flux density of f^, = 2.51 ±0.05 mJy (Figure 1). We subsequently began an intense follow-up 
campaign to study the temporal and spectral evolution of the radio emission. 

Radio data were collected at 1.43, 4.86, 8.46, 15.0, 22.5 and 43.3 GHz spanning 2003 
March through 2005 October and are summarized in Table 1. All VLA observations were 
taken in standard continuum observing mode with a bandwidth of 2 x 50 MHz. At 22.5 
and 43.3 GHz we included referenced pointing scans to correct for the systematic 10-20 arc- 
sec pointing errors of the VLA antennas. We used 3G48, 3G147 and 3G286 (J0137±331, 
J0542±498 and J1331±305) for flux calibration, while phase referencing was performed 
against calibrators J0407-330 and J0453-281. Data were reduced using standard packages 
within the Astronomical Image Processing System (AIPS). No diffuse radio emission from 
the host galaxy was detected in any of our observations. Flux density measurements were 
obtained by htting a Gaussian model to the SN. In addition to the rms noise in each mea¬ 
surement, we include a systematic uncertainty of 2% due to the uncertainty in the absolute 
flux calibration. 

The SN 2003bg radio light-curves are compiled in Figures 2 and 3, spanning ~ 10 to 
~ 1000 days after the explosion. The radio evolution of SN 2003bg is characterized by several 
achromatic short timescale {5t/t < 1) variations. As will be discussed in §8, light-curve 
variations have also been observed for ~ 50% of well-sampled radio supernovae including 
Type HL SN 1979G (Weiler et al. 1991), Type Ic SN 1998bw (Kulkarni et al. 1998; Weiler, 
Panagia & Montes 2001), and more recently for Type Hb SN2001ig (Ryder et al. 2004). 


^The Very Large Array and Very Long Baseline Array are operated by the National Radio Astronomy Ob¬ 
servatory, a facility of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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2.1.1. Radio Spectrum 

Shown in Figure 4 are the spectral indices between each of the adjacent radio frequencies. 
A spectral turnover is clearly observed for the lower frequencies at early time; at later time 
the spectral slope converges to /? ~ —1.1. The observed spectral shape is overall consistent 
with a non-thermal synchrotron spectrum with a signihcant absorption component that 
suppresses the low frequency emission. In the case of radio supernovae, a low frequency 
turn-over may be attributed to internal synchrotron self-absorption (SSA; Chevalier 1998) 
or free-free absorption (FFA; Weiler et al. 1986) from the external medium. The observed 
spectral index of the optically-thick spectrum enables us to distinguish between these two 
processes. For radio supernovae dominated by SSA, the radio spectrum is approximated by 
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( 1 ) 


where Vp is the peak spectral frequency, is the flux density at z/p, and /dthin is the optically- 
thin spectral index (Chevalier 1998). On the other hand, for FFA dominated emission, the 
optically-thick radio spectrum steepens to 
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where Ki and K 2 are normalization constants (Weiler et al. 1986). 

We £t these functions to two well-sampled radio spectra collected on 2003 March 17 
and 2004 February 8 (t ~ 23 and 351 days). As shown in Figure 5, the optically-thick 
spectral component in each epoch is best described with a spectral index of /3 ~ 2.5, and 
thus entirely consistent with the index predicted by the SSA dominated model. Specifically, 
we find the following best-fit SSA parameters: ^110 mJy (50 mJy) and Up ^ 25 GHz 

(2.2 GHz) with Xr ~ 1-5 (4.5) for the early (late) spectrum. The FFA dominated models 
produce steeper optically-thick spectra (/? > 2.5) and thus provide a significantly worse fit 
to the data. We find the following best-fit FFA parameters: Ki ^ 550 (35) and 7^2 ~ 11 
(0.12) with Xr ~ 180 (8-2) for the early (late) spectra. Here we have adopted /dthin ~ —1.1 
for both models, consistent with the observed optically-thin spectral indices from t ~ 10 to 
1000 days (Figure 4). Glearly the SSA model fit is preferred over the FFA fit in both epochs. 
We conclude that internal SSA dominates the absorption on the timescale probed by the 
radio data. 
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2.1.2. Radio Polarization of SN 2003bg 

Since synchrotron emission is inherently polarized, the magnetic held geometry of radio 
SNe can be directly probed through polarization measurements. In order to measure the 
polarization for SN 2003bg, we carried out a full polarization calibration of our 8.46 GHz 
VLA observations taken on 2003 October 7.4 UT {t ~ 227 days). This run was chosen for 
polarization calibration based on the high dynamic range of the radio image. To correct 
for the instrumental polarization, we observed calibrator J0403+260 over a wide range in 
parallactic angle and computed the leakage terms with AIPS task PCAL. We used 3C147 to 
calibrate the absolute polarization angle. 

We do not hnd any signihcant polarization for SN 2003bg. In particular, we place a 
limit on the linear polarization intensity (Jpoi = where Q and I are the Stokes 

parameters) of Jpoi < 0.071 mJy, corresponding to a fractional polarization limit of < 0.8% 
(3 (t). In comparison with other radio SN linear polarization measurements this is one of the 
deepest limits obtained to date: a factor of ~ 4 deeper than the limit for SN 1993J (< 3.3%, 
Bietenholz, Bartel & Rupen 2003) and comparable to that for SN 1979C (< 1%, Weiler et 
al. 1982). 


2.2. X-ray Observations with Chandra 

SN 2003bg was observed with the Chandra ACIS-S detector on 2003 March 24.5 UT 
(t ~ 30 days after the explosion) for 50.5 ksec. Pooley & Lewin (2003) reported the detection 
of an X-ray source at a(J2000) = 04^'10“59.42^ 5(42000) = -31°24'50.3" (±0.5 arcsec in 
each coordinate), coincident with the optical SN position. The source was re-observed with 
ACIS-S on 2003 June 22.3 UT {t ~ 120 days) for 40.5 ksec to search for variability between 
the two epochs. 

We retrieved the SN 2003bg data from the Chandra data archive^ and reduced them 
following the CIAO science threads^. The source is clearly detected in both epochs and we 
note that no diffuse emission is detected from the host galaxy. For spectral extraction of 
the SN emission, we adopt a source aperture of 4.92 arcsec and a large background region 
located ~ 20 arcsec from the source. After subtracting off the background emission, we 
measure count rates of 0.0117 ± 0.0005 and 0.0020 ± 0.0003 cps (0.3-10 keV) in the hrst and 
second epochs, respectively. 


^http://cda.harvard.edu/chaser/ 

^http: / / cxc.harvard.edu/ciao/threads / 
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Pooley & Lewin (2003) report that the first epoch data can be reasonably £t with 
several different spectral models and propose a MEKAL hot plasma model for the best fit. 
Given that the small number of counts prevents the model fits to be distinguished, we fit 
the extracted SN spectra with only two basic models: absorbed power-law and thermal 
bremsstrahlung. Table 2 lists the parameter values for our resulting spectral fits, where we 
have adopted both fixed and variable A^h values for comparison. Schlegel, Finkbeiner & 
Davis (1998) dust maps give E{B — V) = 0.022 mag for the position of SN 2003bg. Using 
the standard conversion of Predehl & Schmitt (1995), we find Ah = 1-3 x 10^° cm“^ for the 
Galactic column density along our the line-of-sight. The relatively few X-ray counts enable 
satisfactory fits ~ 1) for each of the models. 

Given the significantly fewer counts in the second epoch, we fit the spectrum with a 
fixed (Galactic) absorption. Moreover, for the power-law fit we adopt the best-fit spectral 
index from our first epoch. The resulting flux values appear in Table 2 for both power-law 
and thermal bremsstrahlung fits. As in the case of the first epoch, both models provide an 
equally good representation of the data. 

In Figure 6 we show the absorbed power-law model fits for Epochs 1 and 2, assuming a 
fixed Galactic Nh- Using the associated spectral parameters in Table 2, we find unabsorbed 
flux values of Fx ~ (9.3 ± 0.4) x 10“^^ and (1.2 ± 0.4) x 10“^^ erg cm“^ s“^ (0.3-10 keV) 
for these fits, respectively. 

The observed temporal evolution is Fx oc with ax ~ —1.5 ± —0.3 and the implied 
spectral indices between the optically thin radio and X-ray band are nearly constant at 
Prx ~ —0.94 and 0.89, for the first and second Chandra epochs, respectively. At a distance 
of 19.6 Mpc, the observed X-ray luminosity values are thus Lx ~ (4.3 ± 0.2) x 10^® and 
(5.5 ± 1.9) X 10^® erg s“^, placing SN 2003bg among the most X-ray luminous SNe ever 
detected and a factor of ~ 10 fainter than SN 1998bw on a comparable timescale. 


3. Preliminary Constraints 

Readhead (1994) showed that there is an upper limit to the brightness temperature, 
Tg, for SSA dominated radio sources. This limit corresponds to the scenario where the 
fractions of post-shock energy density in relativistic electrons (ce) and magnetic fields (e^) 
are in equipartition, and the shocked electrons are accelerated into a power-law distribution, 
A oc 7 “^, above a minimum Lorentz factor, Xm- Under these assumptions, the equipartition 
brightness temperature, Tep = C{2nk9'^ph'p) Ri 5 x 10^° K, defines the upper bound on 
Tb (Readhead 1994; Kulkarni et al. 1998). Here, is the equipartition size (radius) of the 



radio source. It is noted that Tep depends only weakly on the observed peak frequency and 
flux density. 

Under the assumption that Tb = Tgp, we can estimate the radius and energy of the 
radio emitting material for SN 2003bg. Our first VLA epoch in which the spectral turnover 
is observed is 2003 March 29 {t ~ 35 days). In this epoch we observed a peak flux density of 
fup ~ 85 mJy at Up ~ 22.5 GHz and an optically thin spectral index of P = —(p—1)/2 si —1.1 
(Figure 4). As shown in Figure 5, on the timescale probed by our observations the SN 2003bg 
radio spectrum is well £t with a SSA model and we therefore identify the spectral peak as the 
synchrotron self-absorption frequency, Ua (defined as the frequency where the optical depth 
due to SSA is unity; tssa = !)• Using Equations 2 and 3 of Soderberg et al. (2005a) (hereafter 
S05) the observed source properties i^p, P, and d) imply an equipartition radius of 
6ep ~ 40 /ias (rep ~ 1.2 x 10^® cm) and a total energy for the radio emitting material of 
Eep ~ 1.7 X 10^® erg, where we have assumed Cg = = 0.1. The average velocity of the radio 

shell is therefore roughly Ugp ~ 0.13c and the magnetic field is B^p ~ 2.8e]/^i7ip^ rep^'^^ ~ 0.9 
G. 

A comparison with the ejecta velocities compiled for core-collapse SNe (Ghevalier 1998) 
shows that SN 2003bg resembles more closely a Type Ibc (typical ejecta velocities, v ~ 
0.1c) rather than a Type II event (u ~ 0.01c). Moreover, the inferred ejecta velocity for 
SN 2003bg is a factor of ~ 10 slower than the mildly relativistic speed inferred for SN 1998bw. 
This analysis emphasizes that broad optical absorption lines cannot be used as a proxy for 
relativistic ejecta, consistent with the radio analysis for broad-lined SNe 2002ap (Berger, 
Kulkarni & Ghevalier 2002) and SN2003jd (Soderberg et al. 2006). 

Our preliminary constraint on the minimum energy for SN 2003bg places it among the 
most energetic radio SNe Ibc ever observed, second only to SN 1998bw and comparable to 
SN 2003L (S05). We note that additional absorption processes (e.g. FFA) and departures 
from equipartition (e.g. SN 1993J; Fransson & Bjornsson 1998) increase the energy budget 
further. 

Equipartition analysis may also be used to roughly constrain the characteristic syn¬ 
chrotron frequency, Equating particle kinetic energy across the shock discontinuity, S05 
show that Ura ~ 9.2 x 10^el{v/cpB GHz which gives ~ 0.02 GHz for Cg = 0.1 with the 
values of Ugp and Bep given above. At t ~ 35 days, Um is therefore already below our radio 
observing band and the synchrotron break frequencies are ordered such that Um < i^a- 

Finally, we use these equipartition results at f ~ 35 days to estimate the mass loss rate 
of the progenitor star. From Ghevalier (1998), the mass loss rate of the star, M, can be 
derived from the post-shock energy density in magnetic fields: 
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where we have assumed a stellar wind-blown medium (n oc r ^). We hnd M ^ 6.1 x 
10“^ Mq yr“^ (assuming a Wolf-Rayet wind velocity of = 10^ km s“^). 

For comparison, we estimate the mass loss rate assuming FFA dominates the absorp¬ 
tion. In this scenario, the shock radius is larger than the equipartition estimate and so our 
equipartition parameters provide only a lower limit on the mass loss rate of the progenitor 
star. Using Equation 2.3 of (Fransson, Lundqvist & Chevalier 1996) we hnd: 
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where tffa is the optical depth to FFA processes and T 5 is the temperature normalized to 
10^ K. Here we have adopted the notation lO^Qx = Q where M is given in units of Mq yr“^ 
and Vw is given in km s“h Setting tffa = 1 and T 5 = 1 we hnd > 4.5 x 10^ 

corresponding to a mass loss rate of 0.045 Mq yr“^ {v^ = 10^ km s“^) at t ~ 35 days. This 
is three orders of magnitude larger than that derived from the SSA interpretation. Moreover, 
this would imply that the total energy of the radio emitting material (a proxy for the kinetic 
energy of the fastest ejecta, v > 0.13c) exceeds 1.5 x 10^^ erg. Since hydrodynamic collapse 
distributes the ejecta kinetic energy as Ek oc v~^ (Matzner & McKee 1999), this would 
imply that the optical data (with v ~ 0.06c) have Ex > 10 ^^ erg, larger than any other core¬ 
collapse SN to date. This result, taken together with the fact that a SSA model provides a 
signihcantly better ht to the radio spectra, leads us to conclude that the observed spectral 
turn-over is due to internal synchrotron self-absorption. 


4. Synchrotron Self-Absorption Model 

In modeling the radio light-curves of SN 2003bg, we adopt the formalism of S05 where 
we presented a rigorous formulation of the temporal and spectral evolution of synchrotron 
emission arising from sub-relativistic supernova ejecta. S05 show how the observed SN radio 
emission spectrum at any single epoch is determined by three parameters: Cf,Cr and 
Here, Cf and Cr are normalization constants of the flux density and optical depth with 
cgs units of g s^/^ and respectively. The parameters Cf, Cr and Um are in turn 

determined by the values of four physical parameters: the magnetic held, B, the shock 
radius, r, the minimum electron Lorentz factor, 7 ^, and the ratio ^ With four 
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physical parameters (-B, r, 7 m, 'S) and only three constraints (C/, Cr, r'm,o), we must assume 
an additional constraint to hnd a unique solution. This constraint is obtained by adopting 
a value for 5- By inverting the equations for Cj, Cr and z/m, S05 derive the following 
expressions for 5, r and 7m: 



.0 X 10-« (2 + V/' 

r J / VHz/ 

).3 X (2 + pY^^^p - 2)-Bi7^V17^-i/i 
f Cr /^\-(P-2)/34 
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( 5 ) 

( 6 ) 

( 7 ) 


Here, d is the distance to the supernova and p characterizes the the thickness of the radiating 
electron shell as r/p. The number density of the synchrotron emitting electrons then follows 
directly as 


p-2B^ g 
P - 1 StT nieCpm 


( 8 ) 


where it is assumed that the contribution from electrons in a thermal distribution is negli¬ 
gible. The circumstellar density can be expressed in terms of a progenitor mass loss rate 


StT 2 

M = —Uenipr Vw 

p 


( 9 ) 


by adopting a value for the wind velocity. The total energy of the radio emitting material 
at a given epoch is then given by 


E = 

p ee Stt 


( 10 ) 


where it should be noted that E depends not only on 5^ but also on an assumed value of Cg. 
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4.1. Hydrodynamical Evolution of the Ejecta 

As discussed by Chevalier (1996), several models have been proposed for hydrodynamic 
evolution of the sub-relativistic supernova ejecta. Based on separate assumptions, these 
analytic models constrain the temporal behavior of the shock radius, magnetic held, electron 
Lorentz factor, and energy partition fractions. Here, we adopt the standard SSA model 
(Model 1 of Chevalier 1996) for the hydrodynamical evolution of the ejecta. This model 
assumes that the evolution is self-similar across the shock discontinuity, and thus r oc 
with ar = (n — 3)/(n — s) where n characterizes the density prohle of the outer SN ejecta 
(p oc r“”) and s characterizes the density prohle of the radiating electrons within the shocked 
circumstellar material (ug oc r~^). In addition, the standard SSA model assumes that the 
magnetic energy density {Ub oc B'^) and the relativistic electron energy density {Ue oc rig'-^rn) 
scale as the total post-shock energy density {U oc Uev'^). In this scenario, the magnetic held 
is amplihed by turbulence near the shock discontinuity, implying hxed energy fractions, 
and e^, and thus a constant value of ^ throughout the evolution of the ejecta. 

The evolution of the magnetic held is determined by the CSM density prohle and the 
expansion of the shock radius: B oc with as = oiri^ — s)/2 — 1. We note that for a wind 
stratihed medium, s = 2 and as = —1. The minimum Lorentz factor evolves as oc 
with = 2(ar — !)• From the scalings of Equation 10 it follows that E oc and 

therefore the total post-shock energy increases with time until the swept up circumstellar 
mass is comparable to rest mass of the SN ejecta. 


5. SSA Model Fit for SN 2003bg 

As shown in Figure 2 and briehy discussed in §2, the radio light-curves for SN 2003bg are 
characterized by achromatic short timescale variability, inconsistent with the strict power-law 
evolution prescribed by the standard model. In an ehort to model the overall evolution of the 
radio ejecta we therefore divide the multi-frequency light-curves into four subsections, each 
dehned by the observed time of abrupt variation. We then apply the SSA model described 
in S4 to each subsection. 


5.1. The Synchrotron Self-absorption Peak 

Figure 3 shows that the first peak observed for each of the radio frequencies is chromatic, 
and therefore produced by a cascading spectral break. We dehne the time range for this 
chromatic evolution phase to be t < 110 days, extending until just before the hrst achromatic 
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variation is observed. As discussed in §2 the radio emission is well described by a SSA 
dominated spectrum and we therefore attribute this chromatic subsection to the passage 
of i/a through the observed frequencies. The evolution of the spectral indices (Figure 4) 
indicate that the self-absorption frequency cascades as oc consistent with typical 

values observed for other SSA dominated radio supernovae (Chevalier 1998). 

Using our multi-frequency radio data collected during the time range of the first peak 
and adopting 5 = 1, we £t for the constants Cf and as well as the temporal indices 
ar and As discussed in Section 3, is estimated to be below the radio band during 
our observations and thus is unconstrained by the data. We therefore estimate z/^.o ~ 0.1 
GHz at our chosen reference time of to = 10 days^, consistent with the preliminary value 
calculated in §3. We hnd a best-£t solution {Xr ~ 46; dominated by interstellar scintillation) 
for parameters: p ~ 3.2, C/ ~ 3.1 x 10“^^ g s^/^, Cr ~ 6.2 x 10^® s“^'®, ~ 0.8 and ~ —1. 

These values imply s ~ 2, n ~ 7 and a.y ~ —0.4. As explained within the Appendix of SOS, 
we parameterize the sharpness of the z/^ spectral break with ( = [0,1] and hnd ( ~ 0.6 for 
our hnal solution. This SSA model provides a reasonable £t to the chromatic subsection of 
the light-curves as shown in Figures 2 and 3. 

With this SSA model £t, the physical parameters of the ejecta and GSM are uniquely 
determined by Equations 5 - 10. Chevalier (1982) show that for s ~ 2 and n ~ 7, ?7 ~ 4. 
Adopting this value, we hnd that the expansion of the ejecta is described hj r ^ 6.2 x 
10^^(t/to)°'® cm, and the average velocity is F ~ 0.24(t/to)~°'^c. The radio ejecta therefore 
expand with a modest sub-relativistic velocity comparable to that observed for SN 2002ap 
(Berger, Kulkarni & Chevalier 2002) and signihcantly slower than SN 1998bw (bulk Lorentz 
factor, F ~ 2 on a similar timescale (Kulkarni et al. 1998; Li & Chevalier 1999). Adopting 
Ce = es = 0.1 we hnd that the total energy of the radio emitting material is given by 
E ^ 7.3 X erg, comparable to that of SN 1998bw. The magnetic held evolves 

radially as H ~ 4.9(r/ro)“^'^® G, slightly steeper than that observed for Type Ic SNe 2003L 
(SOS), 2002ap (Berger, Kulkarni & Chevalier 2002) and the Type Ilb SN1993J (Fransson & 
Bjornsson 1998). The electron number density and associated mass loss rate are thus given 
by He ~ 2.2 X 10®(r/ro)~^ cm“^ and M ^ 1.4 x 10“"^ Mq yr“^ where we have assumed a 
Wolf-Rayet wind velocity of = 1000 km s“^ and adopted a nucleon-to-electron density 
ratio of two (appropriate for WR winds). This mass loss rate is consistent with the values 
observed for Galactic Wolf-Rayet stars (Cappa, Goss & van der Hucht 2004). 


■^Here (and throughout) we use subscript “0” to denote the values of parameters at the reference time. 
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5.2. Subsequent Light-curve Variations 

Shortly after the observed SSA peak cascades through 8.46 GHz, an abrupt rise was 
observed at the optically thin frequencies. By t ~ 120 days the radio emission reached an 
achromatic second peak. Subsequent achromatic light-curve variations were later observed 
at t ~ 300 and 600 days and are most pronounced at frequencies above Va- As shown in 
Figure 4 these variations were associated with abrupt changes in the spectral indices, most 
notably for those indices which straddle Va- 

Possible causes for these abrupt light-curve variations include energy injection from a 
central engine, interaction of the reverse shock with density structures in the ejecta, and 
interaction of the forward shock with GSM density enhancements. Kulkarni et al. (1998) 
invoke energy injection to explain the achromatic second peak in the radio light-curves of 
SN 1998bw at t ~ 40 days. This interpretation was supported by the inferred mildly rela¬ 
tivistic ejecta speeds and associated prompt gamma-ray emission, suggestive that SN 1998bw 
was powered by a GRB-like central engine. In the case of SN 2003bg, however, the inferred 
ejecta speed is merely sub-relativistic and therefore there is no evidence for an engine-driven 
explosion. We conclude that the variations observed for SN 2003bg cannot be attributed to 
long-lived engine activity. 

Mioduszewski, Dwarkadas & Ball (2001) show that the late-time flattening observed for 
the SN 1993J radio light-curves can be produced through the interaction of the reverse shock 
with density inhomogeneities within the ejecta. While the 1-D hydrodynamical simulations 
provide a reasonable £t to the smooth evolution of SN 1993J at t > 1000 days, it is not 
clear that they could produce the abrupt rise at t ~ 100 days observed for SN 2003bg. 
Moreover, any abrupt features in the ejecta would likely be smoothed by instabilities in a 
3-D simulation. 

In the third scenario, the interaction between the forward shock and strong density 
variations in the circumstellar material produces abrupt changes in the temporal evolution 
of the radio light-curves. For a wind-stratified medium with constant progenitor mass loss 
rate, the circumstellar density is predicted to follow a smooth profile. Abrupt deviations 
from a wind density profile have been inferred for several radio supernovae (e.g. SN 1979G; 
Weiler et al. 1991) and are argued to be the result of a variable mass loss history and/or the 
effects of a binary companion on the structure of the GSM. These density induced variations 
are more pronounced at optically thin frequencies than for optically thick. Here we show 
that the light-curve variations observed for SN 2003bg are similarly due to GSM density 
enhancements. 

We model the SN 2003bg radio light-curves with a simple density enhancement model 



under the following assumptions: (a) the density variations can be approximated by modest 
abrupt jumps in the CSM such that the density prohle returns to Ue oc after each 
enhancement, (b) the hydrodynamical evolution of the ejecta returns to the self-similar 
solution after each jump, thereby maintaining the same values of cts and both before 
and after each density enhancement, (c) the microphysics of the shock are not affected by 
the density jumps and therefore Cg and eb remain constant throughout the evolution of the 
ejecta, (d) the bulk motion of the ejecta is not significantly affected by the density jumps since 
the slower, more massive material behind the shock continues to plow forward, and (e) the 
minimum Lorentz factor of the radiating electrons is not affected by the CSM enhancements, 
although the number of electrons participating in the radio emission increases. 

Maintaining a self-similar evolution throughout the density variations requires that the 
magnetic held is enhanced by each of the CSM jumps. Moreover, since Ce and eb are as¬ 
sumed to be hxed constants, the additional thermal energy produced by these enhancements 
causes an increase in the total energy of the radio emitting region. Therefore, in htting the 
subsequent peaks in the SN 2003bg radio light-curves, we assume that the evolution of r and 
7 m remain effectively unchanged while B and E vary according to the strength of the CSM 
density enhancements. 

We recognize that this simple density enhancement model is not fully consistent with 
the hydrodynamic evolution of the ejecta. In fact, a small change in the radial expansion 
is expected as the shock wave reaches the higher density region. For the n ~ 7 self-similar 
solution adopted here and the modest density enhancements we infer, we estimate this to 
be minor effect (~ 10%). Therefore, the variations in ejecta parameters {B and E) are 
primarily the result of the CSM density enhancements rather than the minor variations in 
radial expansion. However, the series of density jumps that we infer cause the evolution 
of the radio ejecta to gradually deviate from the initial radial evolution, making our model 
less accurate at very late time (after several years). In addition, we note that the outflow 
does not immediately settle on the new self-similar solution following each density variation, 
but instead settles on a timescale comparable for the radius to double (i.e. the ratio of the 
settling time to the age of the blastwave is of order unity). Thus the sharp features appearing 
in our light-curve hts are merely an artifact of our model assumptions. 

As discussed in §4, Cf, Cr and are functions of just three physical parameters: r, 
7m and B. Of these, only B responds signihcantly to the CSM density variations. From the 
Appendix of SOS we have C/ oc oc and z/m oc B. By fitting for B at each 

density variation, the new values of C/, Cr and z/m are thus directly determined. 

The observed timescales and estimated circumstellar radii for each density enhancement 
are given in Table 3. We note that we do not £t the short rise time of each phase during 
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which the ejecta are settling to the new self-similar solution. As shown in Figure 3, we 
find reasonable fits for each of the light-curve phases by invoking modest (factor of < 2) 
density enhancements. In Table 3 we list the values of Cf, and and the associated 
enhancement factors for the density, magnetic field and ejecta energy for each phase. It is 
interesting to note that the strength of the enhancements decay with time. 


6. Radio Polarization 

As discussed in in §2, our VLA observations constrain the linear polarization of the 
radio emission to less than ~ 0.8%. Although the synchrotron emission from the radio 
ejecta is inherently polarized, the observed polarization level can be significantly suppressed 
for two reasons. First, since the SN 2003bg radiosphere is not resolved in the VLA data our 
polarization limit applies to the integrated polarization of the ejecta emission. As shown 
for supernova remnant Cas A, the integrated polarization can be significantly suppressed 
due to the approximate circular symmetry of the radio ejecta (Rosenberg 1970). Second, 
the radio emission can be depolarized due to internal Faraday dispersion within the ejecta 
(Burn 1966; Chevalier 1982). For the density and magnetic field values we derive from our 
SSA model, we predict the integrated linear polarization of the SN 2003bg radio emission to 
be negligible, consistent with the observations. 


7. Modeling the X-ray emission 

X-ray emission in supernovae can be produced by three processes (Fransson, Lundqvist 
& Chevalier 1996): [1] non-thermal synchrotron emission from radiating electrons, [2] thermal 
(free-free) bremsstrahlung emission from material in the circumstellar shock and/or the ejecta 
reverse shock, and [3] inverse Compton scattering of photospheric photons by relativistic 
electrons. We examine each of these scenarios in the context of the bright X-ray emission 
observed for SN 2003bg. 


7.1. Synchrotron Emission 

As shown in Table 2, the spectral index between the optically thin radio and X-ray 
emission is /? ~ —0.9 for both Chandra epochs. Assuming a constant injection spectral 
index from the radio to X-ray bands we extrapolate the observed radio spectral index of 
1.1 to the X-ray and find that the synchrotron emission under-predicts the observed 



16 


X-ray flux by a factor of ~ 10 in both epochs. This discrepancy grows signihcantly larger 
when the synchrotron cooling break, = I'&'Kmf.ce /is included, beyond which 
the spectrum steepens by A/? = —0.5. Using the magnetic held derived in Section 5, the 
cooling frequency at f ~ 30 (120) days is Vc ~ 6.1 x 10^° (2.3 x 10^^) Hz. Including the 
cooling break and extrapolating from Vc to z/x-ray as fv oc the synchrotron hux falls 

hve orders of magnitude below the observed X-ray emission. We therefore conclude that the 
X-ray emission is not dominated by synchrotron emission. 


7.2. Thermal Bremsstrahlung Emission 

In a scenario dominated by thermal bremsstrahlung processes, the X-ray emission is 
produced by the forward shock plowing into circumstellar material and/or the reverse shock 
heating of the ejecta. In this case, the strength of the X-ray emission depends strongly on 
the density of the CSM and unshocked ejecta. The shocked material then cools by free-free 
emission processes. For a wind stratified density prohle, the summed free-free luminosity 
from the forward and reverse shock is generalized (Chevalier & Fransson 2003; Sutaria et al. 
2003) roughly by 


-^X,FS+RS ~ 8-6 X 10^"^CL^FS+RS^-5'^t«,3^1 ^ ^rg s ^ (11) 

where Cl,fs+rs is a constant such that Cl ,fs+rs = 1 + (n - 3)(n - Af/{An - 8 ). From 
our SSA radio model we hnd ar ~ 0.8 and M _5 14 (^^^,^3 = 1), thus Cf^fs+rs ~ 2.8 and 

Lx:,fs+rs ~ 1.6 X 10^'^ erg s“^. The predicted X-ray luminosity is therefore a factor of ~ 10^ 
fainter than the hrst epoch observation. 

On the other hand, if we assume the X-ray emission must be dominated by free-free 
processes, then this implies a significantly higher mass loss rate of M _5 230 (utu ,3 = 1). 

However, we reiterate that a high density model dominated by external absorption processes 
is inconsistent with the observed radio spectrum (Figure 5). A possible compromise could be 
that the X-ray emission is produced by dense CSM clumps distributed with a very low hlling 
factor (/ << 1) such that the radio emission is optically thin. This scenario was invoked for 
the radio luminous Type Hn SN1986J (Weiler, Panagia & Sramek 1990). However, the dense 
CSM clumps inferred for Type Hn’s also give rise to bright nebular Ha emission lines which 
are not present in the nebular spectra of SN 2003bg. Moreover, it would be difficult to sustain 
the ejecta velocities inferred from our radio modeling in such a a high density environment 
(Chevalier & Fransson, in prep). We conclude that the radio and optical emission do not 
support a scenario in which the observed X-ray flux is dominated by thermal bremsstrahlung 
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emission. 


7.3. Inverse Compton Scattering 

X-ray emission from supernovae can also be produced by inverse Compton scattering 
(ICS) of the relativistic electrons by the optical photons. This process was observed to be 
important in explaining the observed X-ray emission for SN 2002ap (Bjornsson & Fransson 
2004) and SN 2003L (SOS) on timescales comparable to the optical maximum. In this scenario 
the ratio of flux measured in the radio and X-ray bands, respectively is given by 


Tradio ^ Cj? 

Fx Cph 


( 12 ) 


where Ub and Up^ are the energy densities in magnetic fields and up-scattered photons, 
respectively. We note that Equation 12 is relatively insensitive to the spectral energy index 
of the electrons. Bjornsson & Fransson (2004) show that 


Uph ^ 0.4 Lboi,42 j cm ^ (13) 

where td is the time in days, Tboi ,42 is the optical luminosity at t = td normalized to 
10^2 gj,g g-i y jg shock velocity. 

At the time at the first X-ray observation, td ~ 30, the optically thin radio flux 
was Fj-adio ~ 2.2 X 10“^^ erg cm“^ s“^ and the unabsorbed X-ray flux was Fx ~ 9.3 x 
10“^^ erg cm“^ s“^; therefore Fj-adio/Fx ~ 0.24. From our SSA radio modeling (§5) we esti¬ 
mate i? fsi 1.6 G and v ~ 0.19c at t ~ 30 days. Based on the optical peak magnitude (Hamuy 
et ah, private communication) we hnd Fboi ,42 ~ 2.3 at td ~ 30, coincident with the time of 
optical maximum. Therefore, Ub = ~ 0.10 erg cm“^, and Cph ~ 0.028 erg cm“^. 

Thus f/s/Cph ~ 0.28, comparable to the observed ratio for Fj-adio/Fx- We conclude that 
inverse Compton emission can produce the bright X-ray luminosity in the hrst epoch. 

The second epoch, however, cannot be explained through ICS. From our radio model 
at td ~ 120, B K, 0.41 G, u 0.15c and FJ-adio ~ 1-5 x 10“^® erg cm“^ s“F The optical 
light-curve implies Fboi ,42 ~ 0.23 at this time. Therefore we estimate Ub/Upy^ ~ 24 and 
a predicted IC X-ray flux of 6.3 x 10“^^ erg cm“^ s“^ which is a factor of ~ 200 fainter 
than that observed. We conclude that the X-ray flux observed at 120 days is dominated 
by a separate process, possibly due to a cosmic ray dominated SN shock which produces a 
late-time flattening of X-ray light-curves in SNe Ibc (Chevalier & Fransson, in prep.) 



8. Radio Supernovae and CSM Density Variations 


We have shown that the abrupt radio light-curve variations observed for SN 2003bg 
can be attributed to density enhancements in the circumstellar medium. In comparison 
with well-sampled radio supernovae compiled from previous studies^ including our ongoing 
survey of Type Ibc events (Soderberg et al, in prep), we hnd that ~ 50% of radio supernovae 
similarly show evidence for abrupt achromatic variations. At the same time, there are several 
examples of well-studied radio supernovae with smooth light-curve evolutions including Type 
Ilb SN 1993J (van Dyk et al. 1994; Bietenholz, Bartel & Rupen 2003) and Type Ic SN 2003L 
(S05) among others. In Table 4, we compile the light-curve properties for radio supernovae 
with abrupt variability. We emphasize that this compilation contains core-collapse SNe 
from nearly all spectroscopic classihcations, including hydrogen-rich Type IIL, Iln and lib 
in addition to hydrogen-poor Type Ibc. 

Table 4 shows that the short timescale variability observed for radio supernovae can 
be qualitatively divided into three types. The most common effect is an abrupt steepening 
of the radio light-curves as observed for SN 1988Z (Williams et al. 2002). On the other 
hand, an abrupt rise in the optically thin flux has been observed for a couple of events (e.g. 
SN 1987A, Manchester et al. 2002). More intriguingly, episodic (and perhaps periodic) light- 
curve bumps are observed for SNe 1979C (Weiler et al. 1992; Montes et al. 1998, 2000), 
2001ig (Ryder et al. 2004), 1998bw (Kulkarni et al. 1998; Weiler, Panagia & Montes 2001) 
and 2003bg, all of which are characterized by modest flux deviations (factor of < 2). 

These three types of radio variability are generally attributed to dynamical interaction 
of the ejecta with irregularities in the circumstellar medium. Abrupt density variations in the 
media surrounding core-collapse supernovae are known to be produced in a variety of ways, 
including: [1] fluctuations in the progenitor mass loss rate and/or wind velocity. This scenario 
is supported through observations of local supernova remnants and massive stars which 
collectively suggest a complex picture for mass loss during the hnal stage of stellar evolution. 
[2] Termination shocks between stellar winds ejected during the evolution of the progenitor 
star may produce abrupt density enhancements at radii of R> 0.1 pc (Garcia-Segura, hanger 
& Mac Low 1996; Garcia-Segura, Mac Low & hanger 1996; Ghevalier, Li & Fransson 2004; 
Ramirez-Ruiz et al. 2005). Gase 2 is observationally motivated by the circumstellar nebulae 
associated with local Wolf-Rayet stars, which are produced by interacting stellar winds (Ghu 
2002). [3] a clumpy stellar wind perhaps resulting from fragmentation of circumstellar shells 
by Rayleigh-Taylor instabilities (Garcia-Segura, hanger & Mac Low 1996). This scenario 
is motivated by narrow Ha line prohles observed for Type Iln supernovae. In the case of 


® http://rsd-www.nrl.navy.mil/7213/weiler/sne-home.html 
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SN1988Z, Chugai & Danziger (1994) propose a shocked dense wind component at i? ~ 

5 X 10^® cm to explain the low velocity line emission. [4] the effects of a binary companion 
on the circumstellar wind profile. This case is observationally motivated by the binary 
companion for Type Ilb SN 1993J (Mannd et al. 2004) and the dense circumstellar ring 
around SN 1987A at i? ~ 8 x 10^'^ cm (Burrows et al. 1995) commonly attributed to 
the effects of a companion star (Podsiadlowski 1992). In this scenario, the binary system 
parameters determine the circumstellar radii at which the progenitor wind is most affected. 
[5] a close binary progenitor system embedded in a common envelope (Podsiadlowski, Joss 

6 Hsu 1992). Given that the rate of SN progenitors in common envelope binary systems is 
predicted to be low (~ 1% of SNe Ibc; Chugai 1997) this scenario is expected to be rare. As 
an example, a common envelope system has been argued to explain the large inferred density 
enhancement at i? ~ 7 x 10^® cm responsible for strong radio (Stockdale et al. 2004), X-ray 
(Pooley & Lewin 2004) and optical (Ha; Soderberg, Gal-Yam & Kulkarni 2004) emission in 
the unusual Type Ic SN 2001em at t ~ 2.5 years (Chugai & Chevalier 2006). 

To summarize, the strength and location of circumstellar density variations depends on 
the method by which the density profile was disturbed. We note that Cases 3 and 5 are 
generally associated with hydrogen-rich events whereas Cases 1, 2 and 4 are applicable to 
all classes of core-collapse supernovae. 

Finally, in light of this compilation it is interesting to note that an abrupt rise in the 
optically thin radio emission of Type Ibc supernovae can also be attributed to the presence 
of an off-axis gamma-ray burst jet (Paczynski 2001; Waxman 2004). As the jet sweeps 
up circumstellar material and decelerates, it eventually undergoes a dynamical transition to 
sub-relativistic expansion on a timescale of t ~ 1 yr (Frail, Waxman & Kulkarni 2000). Once 
sub-relativistic, the jet spreads sideways and the ejecta rapidly approach spherical symmetry. 
The observational signature of an off-axis CRB is therefore an abrupt optically thin rise as 
the jet sweeps through our line-of-sight. Due to the location of the synchrotron spectral 
peak, this signature is is most easily detected in the radio band on a timescale of t ~ 1 to 
10 years (Soderberg et al. 2006). Yet as we have shown here, many radio supernovae show 
evidence for density induced variability on this same timescale. We therefore emphasize that 
abrupt flux variations in radio supernovae may not be used as a reliable proxy for an off-axis 
CRB jet and/or energy injection from a central engine. 


8.1. A Comparison of SN 2003bg and SN 2001ig 

A comparison of the radio light-curves represented in Table 4 shows that the variations 
are generally unique to each supernova. However, as Figure 8 shows, the separation and 



20 


intensity of the observed variations in SN 2003bg and the Type Ilb SN 200lig are strikingly 
similar. By simply scaling our SN 2003bg light-curve hts by a factor of ~ 0.1 we obtain a 
reasonable match to the SN 2001ig radio evolution. With this scaling we derive the following 
physical parameters for the SN2001ig radio emitting ejecta: r si 4.2 x 10^®(t/to)°'® cm, 
V 0.16c(t/to)~°'^, B ^ 3.8(t/to)~^ G and E ^ 1.3 x 10^®(t/to)°'"^ erg for €e = es = 0.1. The 
electron number density is thus given by Ue ~ 3.0 x 10^(r/ro)“^ cm“^ with successive peaks 
of relative intensity comparable to those observed in SN 2003bg (Table 3) and an associated 
mass loss rate of M ^ 8.6 x 10“^ Mq yr“^ [v^ = 1000 km s“^). We note that this mass loss 
estimate is comparable to that reported by Ryder et al. (2004), but they adopt free-free 
absorption for the dominant absorption process and a progenitor wind velocity of 10 km s“^, 
typical of red supergiant Type II SN progenitors. Our SSA modeling, however, shows the 
ejecta velocity of SN 2001ig is consistent with those of SNe Ibc and the mass loss rate is in 
fact comparable to Galactic Wolf-Rayet stars, similar to the case of SN 2003bg. 

Figure 8 shows that strong light-curve deviations are observed for both SNe on timescales 
of ~ 120 and 300 days with evidence for a weaker feature at t ~ 600 days. Making the 
reasonable approximation that the ejecta evolve as r ~ 5 x 10^®(t/to)°'* cm, these variations 
correspond to density enhancements at radial distances of i? ~ 4 x 10^® and 8 x 10^® cm. 
Assuming a progenitor wind velocity of = 1000 km s“^, we infer ejection timescales of 
~ 12 and 24 years before the SN explosion. Moreover, as shown in the Figure 8 the strength 
of the light-curve variations decreases with time implying that the largest enhancements are 
closest to the explosion site. 

In addition to their radio light-curves, the optical evolution for SNe 2003bg and 2001ig 
are also similar. In both cases there was spectroscopic evidence for a thin layer of hydrogen 
near the time of optical maximum (Phillips et al. 2001; Hamuy, Phillips & Thomas-Osip 
2003); however nebular spectra revealed no hydrogen emission in either event (Filippenko & 
Ghornock 2002; Hamuy, private communication)®. This evolution is characteristic of Type 
lib supernovae (e.g. SN 1993J; Filippenko, Matheson & Barth 1994). However, their inferred 
ejecta velocities resemble more closely those of SNe Ibc and therefore imply compact (Wolf- 
Rayet) progenitors. We conclude that these two events are intermediate between Types 
Hb and Ibc, and show more overall similarities to hydrogen-poor events. Moreover, the 
likeness of their observed radio and optical evolution may indicate similar environments and 
progenitor evolutions. 


®Given the similarity between the two SNe, it is natural to question whether earlier observations of 
SN 200 lig would have revealed a Type Ic spectrum as was observed for SN 2003bg at t ~ 6 days past 
explosion. 
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9. The Progenitors of SNe 2003bg and 2001ig 

The overall resemblance of SN 2003bg and SN 2001ig as outlined in the previous section 
suggests that a similar process is responsible for the observed density variations. Based 
on the location (i? > 5 x 10^® cm) and relative intensity (factor of < 2) of the density 
enhancements, we hnd that Cases 1 and 4 (a variable mass loss rate and the tidal effects 
of a binary companion, respectively) provide the most likely explanations for the inferred 
circumstellar structure. 

Noting the abrupt fluctuations in the SN 2001ig radio light-curves, Ryder et al. (2004) 
recently proposed a binary induced pinwheel nebula for the inferred episodic density en¬ 
hancements. In this scenario, the colliding stellar winds produce a smooth dusty plume 
that spirals in the plane the binary orbit. A highly eccentric orbit is required to produce 
structure at well-dehned intervals and a favorable viewing angle is necessary to reproduce 
the observed flux variations. Here the location and intensity of the pinwheel structure is 
determined by the binary parameters. Ryder, Murrowood & Stathakis (2006) argue further 
support for this model based on the discovery of a late-B - late-F supergiant coincident (rms 
uncertainty of 0.03 arcsec, corresponding to 1.7 pc at d ~ 11.5 Mpc) with the SN 2001ig 
optical position which they interpret as the putative binary companion. However, in cases 
of observed pinwheel nebulae (e.g. Tuthill, Monnier, & Danchi 1999), the dust plumes are 
attributed to the strong wind-wind interaction between a Wolf-Rayet and O- or early-B star 
companion. Taken together with the substantial uncertainty in the position of the putative 
companion, this suggests that a binary induced pinwheel nebula for SN 2001ig is not clearly 
substantiated. 

Adopting a pinwheel model for both SNe 2003bg and 2001ig would imply that the 
geometry of their circumstellar nebulae are nearly identical and thus they share remarkably 
similar binary systems. Moreover, this model would imply that our viewing angle with 
respect to their circumstellar nebulae must also be very similar. While this scenario cannot 
be ruled out, it is clearly unlikely. 

Assuming a single star progenitor model, the mass loss rates derived through our SSA 
modeling for SNe 2003bg and 2001ig imply massive Wolf-Rayet progenitors. Broadband 
observations of Galactic WR stars show that the circumstellar media for these events are 
often signihcantly disturbed. In particular. X-ray observations show evidence for colliding 
stellar winds and dense clumps (e.g. Hillier 2003). If the progenitors of SNe 2003bg and 
2001ig were single WR stars then the observed density variations could be produced by a 
series of enhanced mass loss episodes. To be consistent with the radio data, the mass loss 
ejections must ramp up as the star nears explosion. In this scenario, the observed variation 
timescale and the lack of GSM structure at radii inward of i? ~ 4 x 10^® cm may be explained 
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by a quasi-periodic mechanism which drives enhanced mass loss episodes every ~ 12 years. If 
accompanied by wind velocity variations, these episodes could produce shells of circumstellar 
material, overall consistent with the circularly symmetric density structures we infer for SNe 
2003bg and 2001ig. The observed similarity between SNe 2003bg and 2001ig implies that 
these mass loss episodes may be common among some fraction of hydrogen-poor events. 

These clues suggest that episodic variations in the progenitor mass loss rate and/or wind 
velocity provide a reasonable explanation to the radio light-curve variations observed for both 
SNe 2003bg and 2001ig. Given the predicted timescale for the mass ejections, these two SNe 
provide some of the best evidence for unusual mass loss evolution immediately preceding 
the supernova explosion. While the presence of a binary companion cannot be ruled out 
for either of these hydrogen-poor events, additional studies of short timescale variability in 
radio supernovae will enable discrimination between single and binary progenitor models. 


10. Conclusions 

We report extensive radio and X-ray observations of SN 2003bg. The spectroscopic evo¬ 
lution of this unusual supernova shows a transition from a broad-lined Type Ic to a hydrogen- 
rich Type II and later to a typical hydrogen-poor Type Ibc. This evolution strengthens the 
connection between Type II and Type Ibc (including broad-lined) events. We show that the 
extraordinarily luminous radio emission for a SN Ibc is well described by a self-absorption 
dominated synchrotron spectrum while the observed X-ray emission at f ~ 30 days is ad¬ 
equately fit by inverse Compton scattering of the optical photons off of the synchrotron 
emitting electrons. Our radio model implies a sub-relativistic ejecta velocity, v ~ 0.24c, and 
a size of r 6.2 x 10^® cm at to ~ 10 days. This analysis emphasizes that broad optical 
absorption lines do not imply relativistic ejecta. We hnd that the total energy of the radio 
emitting region evolves as i? ~ 7.3 x erg assuming equipartition of energy be¬ 

tween relativistic electrons and magnetic helds (ce = es = 0.1). The circumstellar density is 
well described by stellar wind prohle, oc r“^, with modest (factor of ~ 2) episodic enhance¬ 
ments which produce abrupt achromatic flux variations. We show that free-free absorption 
does not contribute significantly to the radio spectrum and estimate an average mass loss 
rate of M 3 x 10“^ Mq yr“^, consistent with observed values for local Wolf-Rayet stars. 

Comparison with other events reveals that ~ 50% of radio supernovae show similar short 
timescale light-curve variations which are attributable to circumstellar density irregularities. 
This compilation emphasizes that abrupt radio light-curve variations cannot be used as a 
reliable proxy for an engine-driven explosion, including off-axis gamma-ray bursts. 
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Finally, the radio light-curves and spectroscopic evolution for SN 2003bg are strikingly 
similar to those of SN 2001ig, suggestive of a common progenitor evolution for these two 
events. The overall similarity of SNe 2003bg and 2001ig to radio SNe Ibc is suggestive of a 
compact Wolf-Rayet progenitor model. Based on the relative intensity of the inferred density 
enhancements, we conclude that the progenitors of hydrogen-poor SNe 2003bg and 200lig 
experienced quasi-periodic mass loss episodes just prior to the SN explosion. 
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Table 1. VLA radio flux density measurements of SN 2003bg 


Date 

At 

GHz 

-^!z,4.86 GHz 

Fu,8A6 GHz 

-^i/,15.0 GHz 

Fu,22.5 GHz 

Fu,i3.3 GHz 

Array 

(UT) 

(daysl) 

(mJy) 

(mJy) 

(mJy) 

(mJy) 

(mJy) 

(mJy) 

Conflg 


2003 

Mar 4 

10 

— 

— 

2.51 ±0.07 

— 

— 

— 

D 

2003 

Mar 6 

12 

0.42 ±0.28 

— 

3.86 ±0.10 

— 

30.77 ±0.71 

— 

D 

2003 

Mar 17 

23 

0.55 ±0.20 

2.75 ±0.11 

12.19 ±0.26 

47.03 ± 0.98 

106.30 ±2.16 

— 

D 

2003 

Mar 29 

35 

— 

8.87 ±0.20 

24.72 ± 0.50 

62.11 ± 1.26 

85.39 ± 1.72 

64.07 ± 1.95 

D 

2003 

Apr 11 

48 

— 

14.88 ±0.32 

40.34 ±0.81 

71.51 ± 1.46 

74.58 ± 1.53 

31.75 ± 1.43 

D 

2003 

Apr 21 

58 

— 

22.37 ±0.46 

51.72 ± 1.04 

69.31 ± 1.41 

39.91 ±0.89 

20.56 ±5.87 

D 

2003 

Apr 26 

63 

— 

24.62 ±0.50 

49.64 ± 1.00 

64.75 ± 1.32 

32.48 ±0.71 

15.62 ±2.53 

D 

2003 

May 6 

73 

— 

28.85 ±0.59 

46.20 ±0.93 

41.18 ±0.87 

30.53 ±0.67 

10.18 ±2.55 

TA ^ 
U QO 

2003 

May 18 

85 

— 

31.14 ±0.68 

38.648 ± 0.79 

19.39 ±0.54 

9.63 ± 1.13 

— 

DuA 

2003 

May 24 

91 

— 

32.42 ±0.69 

33.85 ±0.71 

17.34 ±0.83 

10.70 ±2.10 

— 

DnA 

2003 

Jun 17 

115 

— 

39.07 ±0.79 

45.74 ±0.92 

32.14 ±0.68 

21.07 ±0.44 

— 

A 

2003 

Jun 31 

129 

9.27 ±0.30 

47.37 ±0.95 

53.94 ± 1.08 

41.02 ±0.84 

30.19 ± 1.70 

— 

A 

2003 

Jul 4 

132 

8.77 ±0.36 

48.81 ±0.98 

54.27 ± 1.09 

40.92 ±0.83 

— 

— 

A 

2003 

Jul 14 

142 

10.10 ±0.48 

52.69 ± 1.06 

54.83 ± 1.10 

39.24 ±0.81 

— 

— 

A 

2003 

Jul 29 

157 

11.27 ±0.39 

52.81 ± 1.06 

48.43 ±0.97 

30.88 ±0.63 

— 

— 

A 

2003 

Aug 2 

161 

11.38 ±0.44 

50.84 ±1.02 

47.43 ± 0.95 

29.17 ±0.61 

— 

— 

A 

2003 

Aug 22 

181 

12.79 ±0.35 

45.19 ±0.91 

35.76 ±0.72 

21.06 ±0.48 

11.00 ±0.61 

— 

A 

2003 

Sep 11 

201 

14.07 ±0.37 

40.10 ±0.81 

31.35 ±0.63 

16.58 ± 0.44 

6.59 ±0.43 

— 

BnA 

2003 

Sep 24 

214 

17.07 ±0.40 

39.65 ±0.80 

28.67 ±0.58 

16.54 ±0.39 

10.13 ±0.48 

— 

BnA 

2003 

Oct 7 

227 

15.94 ±0.34 

36.45 ±0.73 

27.38 ±0.55 

15.58 ±0.35 

7.84 ±0.47 

— 

BnA 

2003 

Oct 22 

242 

19.31 ±0.49 

36.53 ±0.74 

24.57 ±0.50 

— 

9.38 ±0.23 

— 

B 






Table 1—Continued 


Date 

(UT) 

At 

(daysl ) 

Fy^lAZ GHz 

(mJy) 

Fu,A.86 GHz 

(mJy) 

Fu,8A6 GHz 

(mJy) 

GHz 

(mJy) 

Fu,22.5 GHz 

(mJy) 

Fu,43.3 GHz 

(mJy) 

Array 

Config. 

2003 Nov 4 

255 

19.23 ±0.48 

34.56 ± 0.70 

22.30 ±0.46 

12.71 ±0.35 

7.86 ±0.31 

— 

B 

2003 Nov 15 

266 

23.05 ±0.59 

33.37 ±0.68 

21.67 ±0.44 

11.31 ±0.30 

7.79 ±0.25 

— 

B 

2003 Dec 4 

285 

22.99 ±0.54 

33.12 ±0.67 

21.31 ±0.43 

12.02 ±0.28 

7.74 ±0.19 

— 

B 

2003 Dec 19 

300 

24.13 ±0.54 

33.19 ±0.67 

20.88 ± 0.42 

11.67 ±0.28 

7.76 ±0.19 

— 

B 

2004 Jan 14 

326 

24.50 ± 0.59 

32.90 ±0.66 

20.33 ±0.41 

11.36 ±0.28 

6.84 ±0.19 

— 

B 

2004 Jan 25 

337 

28.30 ±0.61 

31.66 ±0.64 

19.85 ±0.40 

9.08 ±0.26 

5.96 ±0.18 

— 

CnB 

2004 Feb 8 

351 

23.88 ±0.62 

30.90 ±0.63 

18.84 ±0.38 

9.65 ±0.26 

6.72 ±0.16 

3.36 ±0.57 

CnB 

2004 Feb 25 

368 

31.99 ±0.76 

28.73 ±0.58 

17.14 ±0.35 

9.72 ±0.25 

5.72 ±0.17 

— 

C 

2004 Apr 2 

405 

— 

— 

14.61 ±0.30 

— 

7.48 ±0.39 

— 

C 

2004 Apr 7 

410 

33.27 ±0.72 

23.88 ±0.50 

14.49 ±0.31 

— 

— 

— 

C 

2004 Apr 21 

424 

27.93 ±0.63 

22.96 ±0.47 

14.16 ±0.29 

6.41 ±0.25 

2.73 ±0.13 

0±0.39 

C 

2004 Apr 31 

434 

24.82 ±0.57 

20.98 ±0.43 

13.25 ±0.27 

5.75 ±0.21 

3.05 ±0.24 

— 

C 

2004 May 2 

435 

— 

— 

13.08 ±0.27 

— 

— 

0.38 ±0.23 

C 

2004 Jun 29 

493 

26.60 ±0.76 

16.20 ±0.34 

10.04 ±0.21 

6.59 ±0.36 

3.78 ±0.51 

1.31 ±0.38 

D 

2004 Aug 8 

533 

23.65 ± 0.74 

14.60 ±0.33 

8.92 ±0.22 

5.72 ± 0.40 

3.25 ±0.58 

— 

D 

2004 Nov 15 

632 

— 

— 

6.23 ±0.41 

4.47 ± 1.08 

2.74 ±0.57 

— 

A 

2004 Nov 29 

646 

23.70 ±0.51 

10.82 ±0.23 

— 

— 

— 

— 

A 

2005 Jan 24 

702 

23.23 ±0.603 

10.21 ±0.23 

6.18 ±0.14 

3.43 ±0.19 

— 

— 

BnA 

2005 Jan 31 

709 

— 

— 

— 

— 

1.86 ±0.10 

— 

BnA 

2005 Mar 19 

756 

— 

— 

4.62 ±0.16 

2.52 ±0.21 

— 

— 

B 

2005 Apr 11 

779 

23.15 ±0.78 

8.49 ±0.21 

— 

— 

— 

— 

B 



Table 1—Continued 


Date 

(UT) 

At 

(days'! ) 

Fu,1A3 GHz 

(mJy) 

-^!z,4.86 GHz 

(mJy) 

-^!z,8.46 GHz 

(mJy) 

-^!z,15.0 GHz 

(mJy) 

Fu,22.5 GHz 
(mJy) 

Fu,A3.3 GHz 

(mJy) 

Array 

Config. 

2005 May 22 

820 

21.43 ±0.93 

7.56 ±0.21 

3.93 ±0.13 

2.69 ±0.85 

0.85 ±0.23 

— 

B 

2005 Aug 12 

902 

21.71 ± 1.42 

8.23 ±0.20 

4.69 ±0.12 

2.56 ±0.23 

2.22 ±0.35 

1.31 ±0.32 

C 

2005 Oct 27 

978 

21.84 ±0.29 

7.83 ±0.08 

4.48 ±0.07 

— 

— 

— 

DnC 


iDays since explosion assuming an explosion date of 2003 February 22 UT. 



Table 2. Spectral models for Chandra ACIS-S Observations of SN 2003bg 


Date 

(UT) 

At 

(days) 

Model A^h^ 

(10^^ cm“^) 

P^ 

kT^ 

(keV) 

X^/dof 

Flux (0.3-10 keV)t 
(erg cm“^ s“^) 

Prx* 

2003 March 24.5 

30 

Power Law 

o 

o 

1.7 ±0.2 

— 

24/20 

(9.3 ±0.4) X 10"^^ 

-0.94 

— 

— 

— 0.054 ±0.033 

1.9 ±0.3 

— 

24/19 

(9.6 ±0.4) X 10-1^ 

-0.94 

— 

— 

Th. Brem. 

0.013t 

— 

4.3 ±0.9 

26/20 

(7.8 ±0.9) X 10-1^ 

-0.95 

— 

— 

— 

< 0.018 

— 

4.5 ±0.9 

25/19 

(7.7 ±0.9) X 10-1^ 

-0.95 

2003 June 22.3 

120 

Power Law 

0.013t 

1.7 

— 

5/7 

(1.2 ±0.4) X 10-1^ 

-0.89 

— 

— 

Th. Brem. 

0.013t 

— 

4.4 ± 1.6 

4/7 

(1.2 ±0.2) X 10-1^ 

-0.97 


lAll errors represent 90% confidence levels. 

^Nii fixed to the Galactic value (see §2.2). 

* Spectral index between optically thin radio and X-ray flux density. 



Table 3. Radio Light-Curve Peaks for SN 2003bg 


Phase 

At 

(days) 

r 

(xlO^® cm) 

C,* 

(xl0“®^ g s^/^) 

(xlO^S S-3-6) 

u 1 
(GHz) 

rie^ 

(jump) 

(jump) 

(jump) 

1 

<110 

<6.2 

3.1 

6.2 

0.11 

— 

— 

— 

2 

120-250 

6.7-12 

2.7 

23 

0.15 

1.8 

1.3 

1.8 

3 

300-500 

14-21 

2.5 

46 

0.18 

1.4 

1.2 

1.4 

4 

>600 

>24 

2.4 

67 

0.19 

1.2 

1.1 

1.2 


^Values computed at the reference time of to = 10 days. 

^Factor by which the parameter is enhanced with respect to an extrapolation of the parameters from the 
previous phase. 



Table 4. Radio Supernovae with Strong Light-Curve Variations 


SN 

Type 

^peak 

(years) 

^peak 

(z/ GHz) 

Variation 

Type 

Ref.t 

1957D 

II 

< 27 

4.9 

Abrupt Steepening 

1,2,3 

1978K 

Iln 

2.6 

4.9 

Abrupt Steepening 

4 

1979C 

IIL 

1.2 

5.0 

Episodic Bumps 

1,5,6 

1980K 

IIL 

0.4 

5.0 

Abrupt Steepening 

1,7 

1986J 

Iln 

3.8 

4.9 

Abrupt Steepening 

8 

1987A 

II 

0.0093 

0.84 

Abrupt Rise 

9,10 

1988Z 

Iln 

2.5 

4.9 

Abrupt Steepening 

11 

1998bw 

Ic 

0.030 

4.9 

Episodic Bumps 

12,13,14 

2001em 

Ic/II 

< 2.5 

4.9 

Abrupt Rise 

15,16 

2001ig 

Ilb 

0.20 

5.0 

Episodic Bumps 

17 

2003bg 

Ic/II 

0.27 

4.9 

Episodic Bumps 

18 

2004C 

Ic 

< 0.24 

4.9 

Abrupt Rise 

19 

2004cc 

Ic 

0.068 

8.5 

Abrupt Steepening 

19 


^References: 1 - Weiler et al. (1986); 2 - Cowan, Roberts & Branch 
(1994), 3 - Stockdale et al. (2005); 4 - Schlegel et al. (1999); 5 - Weiler 
et al. (1992); 6 - Montes et al. (2000); 7 - Montes et al. (1998); 8 - 
Bietenholz, Bartel & Rupen (2002); 9 - Ball et al. (1995); 10 - Chevalier 
(1998); 11 - Williams et al. (2002); 12 - Kulkarni et al. (1998); 13 - Li & 
Chevalier (1999); 14 - Weiler, Panagia & Montes (2001); 15 - Stockdale 
et al. (2004); 16 - Chugai & Chevalier (2006); 17 - Ryder et al. (2004); 
18 - this paper; 19 - Soderberg et al., in prep. 







34 


MCG-05-10-15 SN2003bg 





t 

10 asec 


Fig. 1.— Composite color {g'r') image of host galaxy MCG -05-10-15 taken with the Palomar 
200-inch telescope with the Large Field Camera. White contours map the radio emission 
from SN 2003bg detected with the Very Large Array at 8.46 GHz on 2005 May 31 (VLA 
B-array conhguration). Contours are linearly spaced and correspond to flux density values 
between 0.5 and 4.0 mJy. The SN was located in a star-forming spiral arm of the host 
galaxy at position a(J2000) = 04^10™59.42’^, 5(J2000) = —31°24'50.3" (±0.1 arcsec in each 
coordinate). We note that no diffuse emission from the host galaxy was detected in any of 
our radio observations. 
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Fig. 2.— Radio light-curves of SN 2003bg. From bottom to top, the measurements were 
taken with the VLA at 1.43 GHz (pink), 4.86 GHz (red), 8.46 GHz (orange), 15.0 GHz 
(green), 22.5 GHz (blue) and 43.3 GHz (purple) between 2003 March and 2005 October 
UT (t ~ 10 — 1000 days after the explosion). The data have been scaled to the 8.46 GHz 
measurements by where /3 ~ —1.1 is the observed optically thin spectral index. An 
abrupt, achromatic light-curve deviation is observed at t ~ 120 days and additional, weaker 
deviations are suggested at t ~ 300 and 600 days. The synchrotron self-absorption model 
described in §5 is overplotted (solid lines). An extrapolation of the early £t (t < 110 days) 
is shown for comparison (dashed lines). 



lux Density (mJy) Flux Density (mJy) Flux Density (mJy) 


1.43 GHz 


8.46 GHz 


22.5 GHz 



43.3 GHz 


100 1000 
Time (days) 


100 1000 
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Fig. 3.— Radio light-curves of SN 2003bg were taken with the VLA at frequencies 1.43, 4.86, 
8.46, 15.0, 22.5 and 43.3 GHz between 2003 March and 2005 October UT. The synchrotron 
self-absorption model hts described in §5 are overplotted (solid line). An extrapolation of 
our SSA model with no density variations is shown for comparison (dashed line). 
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Fig. 4.— Spectral indices for the radio emission of SN 2003bg. The temporal evolution of 
1.43/4.86 GHz, 4.86/8.46 GHz, 8.46/15.0 GHz, 15.0/22.5 GHz and 22.5/43.3 GHz spectral 
indices are shown along with our SSA model fits (solid lines). The abrupt variations in the 
spectral indices are overall consistent with our GSM density jump model (§5). 
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Fig. 5.— Radio spectra for SN 2003bg on 2003 Mar 17 {t 23 days; top) and 2004 Feb 8 UT 
(t ~ 351 days; bottom). Detections are indicated by filled points and upper limits as inverted 
triangles. In the later spectrum, we supplement our Very Large Array data (dark grey points) 
with simultaneous observations from the Giant Meterwave Radio Telescope (light grey points; 
P. Chandra, 2005, PhD thesis, http://www.tifr.res.in/~poonam). The low frequency GMRT 
observations at 0.33, 0.62 and 1.28 GHz confirm that the absorbed spectral index is not 
steeper than (3 ^ 2.5 (Chandra & Ray 2005). Both spectra are best-£t with a synchrotron 
self-absorbed spectrum (thick black lines) as described in §2.1.1. For comparison, we show 
best-fit FFA models (dashed black lines) which are clearly inconsistent with the observed 
spectrum. We conclude that internal SSA dominates the observed absorption. 
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Fig. 6.— Chandra ACIS-S spectra for SN 2003bg at t ~ 30 days (upper panel) and at t ~ 120 
days (bottom panel). Overplotted in both panels is an absorbed power-law fit with hxed 
(Galatic) Ah as described in § 2.2 and Table 2. In §7 we argue that the X-ray emission is 
produced by inverse Compton scattering of the optical photons in the hrst epoch. 
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Fig. 7.— Physical parameters for SN 2003bg based on our synchrotron self-absorption model 
described in §4 and applied in §5. Left column; the temporal evolution from f 10 to 2000 
days is shown for the energy of the emitting region, shock radius and average velocity (top 
to bottom). Right column: the radial prohle of the magnetic held, electron number density 
and mass loss rate are shown from r 5 x 10^^ to 4 x 10^^ cm. An extrapolation of our SSA 
model with no density variations is shown for comparison (dashed line). 
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Fig. 8.— Radio light-curves for SN 2003bg (dark grey circles) and the Type Ilb SN 2001ig 
(light grey circles; Ryder et al. 2004) are compared. While SN 2003bg is overall a factor 
of ~ 10 more luminous than SN 2001ig, the shape of the radio light-curves are strikingly 
similar. Specihcally, the SNe show evidence for light-curve variations on similar timescales, 
perhaps indicative of similar CSM density structures and progenitor system evolution. Our 
radio modeling hts (thick black lines) for SN 2003bg are shown at 8.46 (top panel) and 15.0 
GHz (bottom panel). By dimming the SN 2003bg hts by a factor of ~ 10 we obtain a suitable 
match to the SN 2001ig light-curves. 




